Free electron in a laser field: The nonrelativistic solution 
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Schrodinger equation for a charged particle interacting with the plane wave electromagnetic field 
is solved exactly. The exact analytic solution and the perturbative solution up to second order are 
compared. 
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The exact solution of Dirac's equation for an electron in an external plane wave electromagnetic field was first 
obtained by Volkov [lj in 1935 and it's Green function was derived by Schwinger 0. The state of the charged 
particle, known as Volkov state, has been used extensively to explore numerous quantum phenomena such as Compton 
^vq , scattering, photo-ionization, bremsstrahlung processes and Kapitza-Dirac problem (scattering of an electron by a 
standing light wave), etc. However, an exact analytic solution of Schrodinger equation for a charged particle interacting 
, with a plane wave electromagnetic field has never been obtained. The subject of the interaction of the charged particle 
with radiation is investigated using perturbation theory in almost any quantum mechanics textbooks. However, the 
behavior of a charged particle in the strong radiation fields by high-power laser can not be investigated by perturbation 
theory. To deal with this problem, new approximation methods have been developed 0, 0, H, @, 0, B 0| • For example, 
in dipole approximation, electromagnetic field is assumed to be purely time-dependent for the study of ionization of 
t^J- 1 atoms by intense laser pulses. In this approach, the magnetic field is neglected. As opposes to dipole approximation, 
£>• \ the effect of magnetic field can be taken into account by expanding the potential to first order in the space coordinate. 

■ In this paper, we will solve exactly the Schrodinger equation for an electron in a linearly polarized monochromatic 
plane wave propagating in the z-direction. The vector potential is given by A = {A e ik(z - ct \ 0, 0). The laser pulse 

j. ■ with the elliptic polarization will be discussed later. 

where V.A = 0. To solve this equation, let us transform z- coordinate as s = z — ct. Then the time derivative operator 
transforms as d/dt — ► d/dt — cd/ds. The equation (TT]) becomes 

| : - 4?) = -V 2 * + 2,a e- »* + a 2 * , (2) 
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where a is defined as a = — — - and V 2 = d 2 + d 2 + d 2 We introduce the following ansatz; i.e. we look for the solution 

■ nc v ■ 

• '-j of © in the form 

OS' ®(x,y,s,t)= exp (-i-t - ik' x x - ik'y)^(s) , (3) 
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where e and k' x , k' y are all constants. If we substitute into the equation @, we obtain 

-£ +< ¥S + («' + - 7T> + 2ak > «** + " 2 «"**) < 4 » 

Let us make a complex transformation on the coordinate as u = exp(iks). This is reasonable since the Hamiltonian 

is not Hermitian. Using the transformation for the derivative operator — = iku—, the equation (|U) becomes 

as du 
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mc K y + K x — p-e 

where a = ——, 5 = — — — . To solve (15]), the last transformation is introduced as 

hk k 2 



$(u) =exp(~i— ) v? L{u) , (6) 
k 



2 



where 7 is a constant to be determined later and L{u) is a function depending on u. If 7 satisfies a polynomial 
equation j 2 — 2(77 + <5 = 0, then L{u) satisfies 



"A? + ((27 + 1 - 2a)u - 2i-u>) _ + ^ _ ,_ ( 2 7 + ! _ 2(T ) j ul=0. (7) 

The solution is given in terms of the associated Laguerre functions L 21 2a ,, ( —, — | . Note that a generaliza- 

<r-7-iV-i V * y 

tion of the associated Laguerre polynomial to not necessarily an integer is called the associated Laguerre function. 
Transforming backwards yields the exact solution 

* = exp (-tit - ik' x x - ik' y y - i^ e tk(z - ct) + ijk{z - ct)) L 2 ^ 2 " k , ± Ct) ^ . ( 8 ) 

To determine e and 7, let us set the charge to zero (a — 0). In this case, the free particle solution should be recovered. 
Note that Laguerre function is a number at zero and the constants can be omitted since the free particle solution isn't 
square integrable. 

To obtain the free particle solution, let us choose 7 = —k' z /k. Then, by solving the polynomial equation for 7, we get 
—5 = -r(-r + 2a). Using the definition of 5, we obtain e — hk 'c — - — Vk' 2 + k' 2 + k' 2 ). Defining a new constant 



k' k' h 2 
-r(-r + 2a). Using the definition of S, we obtain e — Tik'c = — , 
k k 2m 
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as E = e + jTikc = - — k' , the free particle solution is recovered. Physically, E is the energy for the free particle. 
2m 



Substituting j,5,e,E in ©, we get 

*( W ,i)= ^ V [-i^t-ik'.T-i^-^\ L- 2 J-*° l^r 1 ) • (9) 
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This is the exact solution of Schrodinger equation with the plane wave electromagnetic field polarized in the x- 
direction. If we set a to zero, then it is reduced to the free particle solution. The first two terms in the exponential 
is the exact solution in the absence of the external field, the other two terms is due to the interaction with the field. 
The solution has been obtained for the linearly polarized electromagnetic field. If the field is elliptically polarized, 
then the term k' x in Laguerre function should be replaced by the linear combination of k x ,k y . 

Finally, a few remarks about the solution ^ is in order. Some approximations can be made for the associated 

Laguerre function L k , . The order of a is only comparable to the factor 1/2 when k « 10 11 m _1 (x-rays). 

Otherwise, one of them can be neglected. Furthermore, if k! is small compared to k, then this term can be neglected, 
too. 

As a special case, perturbative and our non-perturbative solutions should yield the same results. To check the validity 

of our result, let us compare the exact solution ^ with the the one obtained by perturbation theory. Using the series 

r jt c rmi \ T(m + n + 1) ( nx (l-n)nx 2 \ 

expansion ol associated Laguerre function L„ (x) — — — — — 1 — — — — — ... , where 

1 & nK 1 T(m + l)T(n + 1) V m+1 2(m+l)(m + 2) /' 

r is the gamma function, the wave function ([9]) can be expanded in series up to order a 3 as 

* = e -i(K*+Kv) {e- ik >-^ + a ^-HK-V'-i^P* + a * c ^iK-^V-^^ + (10) 



where c\ — — ; 7^7 , C2 = - , 2 -, — z - N i ^, Unl , — , /n 777 ■ The constant in front of this expansion 



2k' x _ (2kk' z + 4k' x A + k 2 (2a-l)) 

k{2k' z + k{2a-l)Y ° 2 ~ 4k 2 (k(a - 1) + k' z )(2k' z + k[2a - 1)) 
can be dropped from the solution since it is not physical. 

Having obtained the series expansion of the exact solution (9), let us review the perturbative expansion for the 
problem. Assuming a is small, perturbation theory states that can be approximated as 

* = -0o + a &iV>i + a 2 b 2 ip 2 + 0(a 3 ) , (11) 

where tf>o is the free particle solution. Substituting the equation pip into the corresponding Schrodinger equation ((T|) 

2 

yields a b\ ( ipx — ihipi) = — Ap ipo for the first order of a. Furthermore, for the second order a 2 , we are left 

2m mc 

2 2 

with the following equation a 2 ^2 (7; — i>i — ihtp2) = b±a — Ap ib\ — — 7 A ipo- By solving these two equations, we 

2m mc 2mc l 
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see that b\ = c\ and 62 = ci- Then, it is concluded that our non-perturbative exact solution is in agrement with the 
solution obtained by using the perturbation theory. 

Perturbation theory fails when applying to the system with strong electromagnetic fields. In contrast, the solution 
([5]) is non-perturbative exact solution and it can be applied without restricting ourselves to the weak laser field. 
The non-perturbative exact solution obtained here can be used in many branches of physics. For example, the recent 
development of high power laser delivering pulses of intensity up to 10 22 W/cm 2 enables us to investigate deeply 
the properties of atoms, molecules, plasmas and condensed matter interacting with super intense laser pulses. With 
rising laser intensity, non-perturbative treatments are required to study the electronic motion of free and bound 
particles in intense laser fields. Analytical studies of strong laser atom interaction have seen much progress over the 
years 1(J 11, 12, 3 EH- A successful theory so called strong field approximation (SFA) was developed by Keldysh, 



Faisal and Reiss [10(. In SFA, the final state of the ionized electron is represented by Volkov solution which describes 
a free electron interacting with laser field. This is because laser field dominates Coulomb interaction. With this 
approximation, the physical phenomenons such as above-threshold ionization and high harmonic generation have 
been explained successfully. Exact solution for Coulomb field is well-known and given in terms of associated Laguerre 
polynomials. The exact solution for an electron interacting with radiation is given in the equation ((9]). Thus, for the 
problem of atom-laser interaction, either field dominates the other can be treated perturbatively. 
The dipole approximation is often be employed in atomic physics. If the Bohr radius ao is very smaller than the 
wavelength of the radiation ao << A, then the phase can be approximated as uit — kz ss uit. The consequence is that 
the magnetic field is neglected, since magnetic field is given by B = V x A. If the laser field is sufficiently intense, 
then a fully relativistic treatment for the dynamics of the electrons is required. Since the wave length is not small, 
the dipole approximation fails for the intense laser fields. 



As it was discussed in |15j . there is, however, an intensity region intermediate between the nonrelativistic and 



relativistic domain. The force on a charged particle due to the magnetic field is related to v/c. This can be seen 
from the Lorentz force formula F = e(E + v/c x B). True relativistic effects, however, are v 2 /c 2 effects. Hence, there 
is an intermediate region where the dipole approximation is no longer valid and the relativistic treatment is not 
necessary. In this region, the particle is described by Schrodinger equation which includes all relativistic effects of 
order v/c. Some methods have been developed to understand the dynamics of laser-matter interaction in this region 
Til 13) E3) 18, The exact solution obtained here can also be used for the domain where v/c but not v 2 /c 2 are 



necessary. 

Finally, we think that the analysis of high harmonic generation 2^, 21 1 can be performed more accurately by using 
the exact solution (O. 
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